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Laterally periodic nanostructures were investigated with grazing incidence small angle X-ray scat¬ 
tering. To support an improved reconstruction of nanostructured surface geometries, we investigated 
the origin of the contributions to the diffuse scattering pattern which is correlated to the surface 
roughness. Resonant diffuse scattering leads to a palm-like structure of intensity sheets. Dynamic 
scattering generates the so-called Yoneda band caused by a resonant scatter enhancement at the 
critical angle of total reflection and higher-order Yoneda bands originating from a subsequent diffrac¬ 
tion of the Yoneda enhanced scattering at the grating. Our explanations are supported by modelling 
using a solver for the time-harmonic Maxwell’s equations based on the finite-element method. 


I. INTRODUCTION 


X-ray scattering is established as a common tech¬ 
nique in nano science. Among the diverse nano-sized 
object spresently under investigation in natural and life 
science^®, the continuously shrinking patterns in the 
semiconductor industry are at the forefront regarding 
their requirements for size reproducibility and regular- 
The latter also differ from many other applications 
which benefit from the single particle diffraction measure¬ 
ments now feasible with free-electron laser sourceJ^^^^. 
Measurements on the length scale of several nm are basi¬ 
cally challenged by the atomic granularity of matter and 
by structures which are not perfect and sharply defined 
anymore. Due to the large photon beam footprint as 
compared to the nm pattern size, grazing incidence small 
angle X-ray scattering (GISAXS)!^ directly yields statis¬ 
tical information for a large structured area on fluctua¬ 
tions like structure roughness. It is a powerful tool for the 
analysis of the morphology and dis tribut ion of nanoparti¬ 
cles on surfaces or buried particle^^. As compared to 
SAXS in transmission geometr^ b^F^ I^ GISAXS also allows 
to study structured surfaces of thick, non-homogeneous 
substrates. 

Recent GISAXS studies on nanoscale line gratings 
mainly focus on the grating diffraction and deliver 
no explanation for the diffuse scattering background. 
Gorresponding theoretical calculations for small angle 
X-ray scattering in the past were mostly don e with t he 
distorted wave Born approximation 
including an analytic expression of the line shape. Arbi¬ 
trarily shaped structures must be discretized and require 
numerical solverJ^. For periodic structures, e.g. grat¬ 
ings, it is well established to model the light scattering 
by numerically solving the time-harmonic Maxwell’s 
equa tions with a higher-order finite-element method 
If file periodic structures are invariant in one 
dimension (along the grating lines), the computational 
domain can be reduced to a two dimensional problem 
which decreases the computational effort significantly 
and allows to calculate also rather large domains as 


compared to the incident wavelength. This helps to 
apply the method to short-wavelengths like in X-ray 
scattering. The computation yields a near field solution 
for the electric field in the structured volume. This can 
be post-processed to obtain the far field solution, i.e. the 
scattering intensities, or other near field parameters like 
the electric field strength along a certain path, e.g. the 
vacuum-material interface. 

The scope of this investigation is to determine the kine¬ 
matic and dynamic scattering processes which cause the 
diffuse background and relate them to geometrical pa- 



FIG. 1. (a): SEM image of a silicon grating sample similar to 
grating No. 5 (see table in Fig.[^. (b): Finite-element model 
employed for the simulations of above grating sample, (c): 
Scattering of a collimated X-ray beam under grazing incidence 
from a lamellar grating oriented parallel ((y^ = 0°) with respect 
to the incident beam. 
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rameters of the scattering surface in order to obtain more 
detailed information on nano-structured surfac e geome¬ 
tries. Yoneda enhanced surface scatter in^^^^^ and sub¬ 
sequent diffraction at a periodic nanostructure give rise 
to the formation of higher-order Yoneda bands in the dif¬ 
fuse scattering background. Meanwhile, the origin of the 
Yoneda bands and additional waveguide effects can be 
explained with a grating effective layer. 


II. EXPERIMENTAL DETAILS 

We investigated lamellar silicon gratings as a proto¬ 
type specimen for state-of-the-art integrated electronic 
elements. The grating structures were prepared by elec¬ 
tron beam lithography and reactive ion etching on a sil¬ 
icon wafer. The grating area measures 1 mm by 15 mm 
with the lines oriented parallel to the long edge. This 
accommodates the elongated footprint of the beam in 
GISAXS geometry. The five samples discussed here are 
listed in the table below Fig. 

The experiments were conducted at the PTB’s four 
crystal monochromator beamlin^^ at the electron stor¬ 
age ring BESSY II, which is equipped with an in-vacuum 
Pilatus IM area detect oi!^. The GISAXS scheme is illus¬ 
trated in Fig. [^: A monochromatic X-ray beam with a 
wave vector ki impinges on the sample surface at a graz¬ 
ing incidence angle ai and an azimutal angle Oi = 0°. 
The elastically scattered wavevector kf propagates along 
the exit angle af and the azimuthal angle Of. 


III. DIFFUSE SCATTER CONTRIBUTIONS 
FROM STRUCTURED SURFACES 

A. Yoneda band - Waveguide effects 

The regular structure of a grating leads to discrete 
diffraction orders on a cone with an opening angle equal 
to the incidence angle ai. Roughness and imperfections 
of the grating give rise to satellite diffraction orders, 
which result from the fabrication process and lead to dif¬ 
ferent diffuse scattering contributions. If the exit angle 
af is equal to the critical angle of total reflection o^c, the 
diffuse scattering from a rough surface is resonantly en¬ 
hanced. This resonant e nhanc ement is often called the 
Yoneda or Vineyard pealP^^ which appears here as a 
bright horizontal line of enhanced diffuse scattering in¬ 
tensity. The mechanism behind this enhancement is vi¬ 
sualized in Fig.j^by showing the electric field intensity at 
the surface for incident angles around the critical angle 
(a). The field intensity, and thus the potential scatter¬ 
ing from the surface, reaches its maximum at the critical 
angle. 

If the surface is capped with a thin layer (b), the 
Yoneda line splits up into a band with pronounced oscilla¬ 
tions between the critical angles of the substrate and the 
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FIG. 2. Calculated electric field intensities at the vacuum 
interface for incident angles ai around the critical angle of 
silicon, (a): The dashed black line represents the surface field 
enhancement at the critical angle for a flat Si-substrate, (b): 
The substrate was capped with an additional layer with re¬ 
duced mass density (42%) and a thickness of 100 nm. (c): The 
field intensity for a grating with rectangular lines at oblique 
azimuthal incidence (cp = 2°). 



c) Measured Yoneda band _ d) Measured Yoneda band 
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Grating 

No. 

nominal 

height 

(nm) 

nominal 

pitch 

(nm) 

nominal 

structure 

density 

sim. 

Qi (deg) 

exp. 

a/ (deg) 

I 

100 

100 

0.35 

0.II7 

0.115(12) 

2 

100 

250 

0.40 

0.125 

0.122(10) 

3 

100 

100 

0.45 

0.135 

0.127(14) 

4 

50 

50 

0.50 

0.147 

0.147(13) 

5 

100 

100 

0.55 

0.144 

0.141(11) 


FIG. 3. The electric field intensity at the vacuum interface 
normalized to the incident plane wave calculated for ai below 
the critical angle of silicon for a rectangular grating with a line 
height of 50 nm (a) and 100 nm (b) as a function of structure 
density. The angular scale is normalized to ac- Blue circles 
in both figures mark the experimentally determined position 
a/ of the lower Yoneda band boundary. The numerical values 
for the simulated angles ai and measured a / are summarized 
in the table. (c,d): Measured Yoneda band characteristic 
obtained from grating No. 4 and grating No. 3. 
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log. intensity (arb. units) 

IQi 10^^ 10^ 10^^ 10^ 



FIG. 4. GISAXS pattern from grating No. 4 obtained at 
9.5 keV and ai = 0.5°, (f = 0° with logarithmic gray scale. 
The range of the scale was adjusted to optimize the visibility 
of the diffuse scattering background. 


tained by matching the calculated and measured (Fig. 
blue circles) lower band boundaries. The sharpness of the 
Yoneda bands for the individual samples varied, resulting 
in different error margins as shown in Fig. 


B. Higher-order Yoneda scattering 

The Yoneda scattering process gives rise to an addi¬ 
tional structure in the diffuse scattering background. In 
contrast to a homogeneous layer, the grating causes also 
higher ordering of the Yoneda band, appearing as bent 
sharp lines originating from the sample horizon and cross¬ 
ing the Yoneda band (Fig. |^. In a first step the Yoneda 
effect enhances the diffuse scattering intensity at exit an¬ 
gles af close to the critical angle and for arbitrary az¬ 
imuthal exit angles. This is the 0^^ order Yoneda band. 
For a periodic structure like a grating, also higher or¬ 
dering occurs if the Yoneda-scattered beam is diffracted 
in a second step at the grating. Clear evidence for the 
role of the grating diffraction in the higher order Yoneda 
bands is the fact that an azimuthal rotation of the grating 
leads to a corresponding shift of these bands preserving 


capping layer caused by the formation of an X-ray stand¬ 
ing wave inside the capping layeJ^. The incident wave is 
coupled to these waveguide modeP^. This effect has been 
used for dep th profiling of thin films or buried nanostruc- 
turefl^E^^. Here, we observed a similar Yoneda band 
in the GISAXS pattern from laterally structured surfaces 
instead of a simple layer system. Electromagnetic field 
computations for a rectangular grating (c) and oblique 
incident angles ai demonstrate that the local field en¬ 
hancement above the grating behaves almost identical to 
that of a top layer with reduced density. For diffra ction 
efficiency measurements this behavior is well knowii^^^ 
from scattering geometries perpendicular to the grating 
lines {ip = 90°). 

Here, we investigated the Yoneda band characteristics 
for several lamellar silicon gratings with different geomet¬ 
rical dimensions (see Table in Fig.[^ close to the parallel 
orientation of the lines with respect to the incoming beam 
{ip = 0°). The number of waveguide modes depends on 
the structure height and density, as illustrated in Fig. 
The electric field at the vacuum interface, following the 
line profile into the grating grooves, was computed for an 
azimuthally rotated rectangular grating {ip = 1°) and in¬ 
cident angles ai around the critical angle of silicon. For 
the grating with 50 nm nominal height and a 1:1 duty 
cycle, measurement and electric field computation con¬ 
sistently yield a splitting of the Yoneda band, illustrated 
in Fig.[^. The gratings with 100 nm nominal height show 
pronounced oscillations inside the Yoneda band very sim¬ 
ilar to waveguide modes (Fig. If)- Variation of the grat¬ 
ing linewidths reveals the high sensitivity of the Yoneda 
band for structure density, i.e. the line-to-pitch ratio, as 
well as height variations. Vice versa, experimental values 
for the linewidth with respect to the pitch can be ob- 


log. intensity (arb. units) 

10^ 10^-5 10^ lO^-^ 10^ 




FIG. 5. (a): The central area under an azimuthal rotation 
(F = 8 keV, ai = 0.8°, p = 0.5°) of the grating No. 4. The 
diffraction cone and higher order Yoneda bands follow this 
rotation, in contrast to the RDS sheets. Red dots illustrate 
the expected reciprocal space locations for the higher orders 
of the Yoneda band, (b): Magnified view of figure (a) around 
the higher-order Yoneda bands. 
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FIG. 6. Magnified view of the diffuse scattering pattern obtained from (a) grating No. 5 = 6 keV, ai = 0.1°, (/? = 0°) and 

(c) grating No. 4 (F^ = 5 keV, ai = 0.15°, = 0°) around the higher-order Yoneda bands. (b,d,e): Computation of the electric 

field intensities for both gratings (see text). The insets show the rectangular grating model including corner rounding (b,d) 
and a trapezoidal model with 80° sidewall angle (e). 


the scattering direction with respect to the grating lines 
(see Fig. [^). We calculated the reciprocal space maps 
for radiation first scattered at af = ac, the critical angle 
of either the substrate or the effective reduced density 
surface layer, and a random azimutal angle Of and sub¬ 
sequently diffracted at the mating structure, following 
the formalism of Yan et for a grating azimuthally 
rotated. The resulting reciprocal space locations for 81 
equally spaced values of Oi, i.e. of the first Yoneda 
scattering event, in the range ±0.4° agree well with the 
observed scattering pattern (Fig. red dots). 

An azimuthal rotation of the grating allows to separate 
the higher-order Yoneda bands from other scatter contri¬ 
butions close to the specular axis. Closer inspection of 
Fig. 1^, where the grating was azimuthally rotated with 
respect to the beam ((^ = 0.5°), reveals a complex inten¬ 
sity distribution of the higher-order Yoneda bands. The 
resonant Yoneda enhancement of diffuse scatter contribu¬ 
tions interferes with the grating structure at the critical 
angles of substrate and effective layer and produce a char¬ 
acteristic scattering pattern which is directly related to 
the geometrical shape of structure. This relation is also 
demonstrated in Fig. with a comparison of the diffuse 
scattering patterns obtained in conical scattering geom¬ 
etry {(p = 0°) from grating No. 4 and 5. The distinct 
differences in the geometrical dimension of the gratings 
produce a significant variation in the diffuse scattering 


patterns. To emphasize the higher-order Yoneda bands 
the gratings were measured in total reflectance with an 
incidence angle ai well below the critical angle of sub¬ 
strate and grating effective layer. 

The extremely large computational domain required to 
represent roughness in a sensible manner renders a direct 
modeling of diffuse scatter with a finite-element method 
impractical. We follow here the same reasoning as for 
the 0^^ order Yoneda enhancement, i.e. we suppose that 
the scatter is intense if the electrical field at the grating 
surface is high. We thus computed the electric fields at 
the vacuum interface of the grating lines for various inci¬ 
dent angles (ai^Oi) for both gratings. In analogy to the 
Yoneda enhancement, the scatter intensity at the corre¬ 
sponding exit angles (af^Of) should follow this distribu¬ 
tion. The computed pattern for a finite-element model 
similar to the SEM cross section image of grating No. 5 
with pronounced corner roundings inside the grooves as 
shown in Fig. is indeed very close to the measured 
scatter distribution in the higher-order Yoneda bands in 
Fig.j^. The simulation of grating No. 4 (cf. Fig. and 
Fig.^) reflects also nicely the measured scatter distribu¬ 
tion. To demonstrate the sensitivity for model variations 
we introduce a reduced sidewall angle (80°) in Fig. 
to obtain a trapezoidal grating line shape. The clearly 
visible difference (cf. Fig. and FigJ^) offers an oppor¬ 
tunity for reconstruction of nm sized surface structures 
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based on the measurements of the diffuse scattering pat¬ 
tern. 


C. Resonant diffuse scattering 

Apart from Yoneda resonant scattering, a palm-like 
structure of stacked sheets around the specular axis (see 
Fig.0 is observed. Simil ar pa tterns can be seen in the 
figures of previous studies^^^ but they were never dis¬ 
cussed nor explained. In contrast to the higher-order 
Yoneda band, the position of these sheets is invariant 
with respect to ai and variations. We attribute these 
sheets to resonant diffu se scattering (RDS) known from 
multilayer system^^^^. There, their periodicity in the 
direction of momentum transfer perpendicular to the sur¬ 
face, represents the layer thickness. For the line 
gratings we observe similar behavior with respect to the 
line height, but with an additional dependence on the 
qy momentum transfer, which bends the sheets. Mea¬ 
surements of the diffraction efficiencies for various inci¬ 
dent angles show an identical modulation for the differ¬ 
ent diffraction orders. Computations of diffraction effi¬ 
ciencies and^Cray diffraction experiments for trapezoidal 
grating^ ^^ l ^^ l ^^ l predict a linear shift of the modulation 
along qz with whose slope is directly related to the 
sidewall angle of the grating structure and valid for trape¬ 
zoidal gratingJ^ni, SEM cross-section images (cf. Fig. Im) 
reveal that the line shape of the gratings investigated 
here deviates from a perfect rectangular or trapezoidal 
profile, which explains the palm-like bending of the RDS 
sheets. 


IV. CONCLUSIONS 

In conclusion, we have presented a full explanation 
of the diffuse scatter contributions and indicated ap¬ 
proaches for numerical modeling which pave the way 
to exploit these effects for reconstruction of surface 
structures. In equivalence to X-ray standing wave 
effects, known from vertically layered systems, a Yoneda 
band including waveguide modes forms at a laterally 
structured surface. Here, the grating structure acts as a 
layer with reduced density and the resulting waveguide 
modes provide a direct measure of the structure average 
density, respectively linewidths for rectangular lines at 
fixed pitch. Higher order Yoneda bands result from 
enhanced surface scattering and subsequent diffraction 
at the grating. Their reciprocal space representation is 
predicted with an analytic approach and a numerical 
computation of the electric field at the surface of the 
nano structure correlating the Yoneda resonant enhance¬ 
ment with the surface figure. The palm-like bending of 
the resonant diffuse scatter sheets is also linked to the 
surface figure, in our example the shape of the grating 
lines. 
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